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, Abstract 

m ■ 

i Supersymmetric QCD corrections to the gluonic production and decay rate of a CP-even 

Higgs boson are evaluated at next-to-leading order. To this aim, we derive an effective 
Lagrangian for the gluon-Higgs coupling. We show that a consistent calculation requires the 
inclusion of gluino effects, in contrast to what has been done previously. The supersymmetric 



Oh! corrections to the giuon-Higgs coupling lead to a modification of the next-to-leading order 

^ I K-factor for the Higgs production rate at the LHC by less than 5%. 

> 

^ ■ 1 Introduction 

It is among the main goals of the CERN Large Hadron Collider (LHC) to observe the Higgs boson 
and study its properties. The most important production mechanisms for a Standard Model 
(SM) Higgs boson have been studied in great detail and are known with fairly high theoretical 
precision (for a review, see Refs. [Bill [3]). The dominant mode is gluon fusion, where, at leading 
order in perturbation theory, the gluons couple to the Higgs boson preferably via a top quark 
loop. A light Higgs boson, Mjj < 130 GeV, would be observed through its subsequent decay 
into photon pairs. For larger values of the Higgs boson mass, H ZZ^*^ and H WW^*^ 
allow for a clear identification of the Higgs signal. 

The next-to-leading order (NLO) corrections to the gluon fusion process in the Standard Model 
are known exactly [H [3 [2] . They were shown to be approximated extremely well by an earlier 
calculation which was based on the heavy-top limit [6l [7j . This is true for a large range of the 
Higgs mass, including even Mh > 2Mt, provided that the mass dependence at leading order in 
as is kept exactly. Recently, the next-to-next-to-leading order (NNLO) corrections in the heavy- 
top limit became available [H [9] (see also Ref. [lOj). They lead to a significant stabilization of 
the theoretical prediction, with a remaining scale uncertainty of around 15%. The SM Higgs 
production rate at the LHC is thus under good theoretical control. 
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In the Minimal Supersymmetric Standard Model (MSSM), several aspects have to be taken into 
account when considering the production of a neutral scalar Higgs boson: For example, for 
large tan/3, the bottom quark can contribute significantly to the gluon-Higgs coupling. This 
effect has been studied at NLO in Refs. [H [HI [5l [2]. A large value of tan/3 also brings in 
the associated production of a Higgs boson with bottom quarks as one of the main production 
mechanisms [HI [131 [H] • 

The main subject of the current paper are the contributions of MSSM particles to the gluon- 
Higgs coupling. At sufficiently small values of tan/3, the dominant effects are due to top-squarks 
and gluinos. The NLO squark contributions have been addressed in Ref. [15j. However, as we 
will show, these results do not correspond to a consistent supersymmetric framework. In fact, 
it turns out that the squark contributions can not be separated from the gluino effects. This 
will significantly affect the final result, most evidently in the appearance of non-decoupling 
logarithmic contributions from the gluino mass. 

Motivated by its success in the SM case, we construct an effective Lagrangian at NLO for the 
gluon-Higgs interaction in the MSSM. This corresponds to integrating out the top quark and all 
supersymmetric particles. Within the framework of this effective Lagrangian, all calculations 
can be performed in standard five-fiavor QCD. They can thus be taken over from the SM case. 

In this paper, we restrict ourselves to the case without squark mixing, and focus on a principle 
study of the method. A more general analysis will be presented in a forthcoming paper. 



2 The effective Lagrangian in supersymmetric QCD 
2.1 General definition 

It has been shown that the radiative corrections to the production rate of a SM Higgs boson 
are well described by an effective theory, where formally the top quark is assumed much heavier 
than the Higgs boson. In practice, it turns out that this approximation works very well even 
way beyond the formal limit Mt > Mh- We therefore assume that a similar behavior will be 
observed in the MSSM, and construct an effective Lagrangian where the top quark, the squarks, 
and the gluino are integrated out. The Higgs boson itself enters only as external particle and 
thus the Lagrangian describing the effective (j)gg coupling at leading order in has the form 

Cy,.^ = -^CW, 0? = iG°.G°'^^ (1) 

where Ci is the coefficient function containing the remnant contribution of the heavy particles 
and Oi is the effective operator, (f) denotes either of the two CP-even Higgs bosons of the 
MSSM. is the gluonic field strength tensor of standard, five-flavor QCD. The superscript 
"0" indicates bare quantities as opposed to renormalized ones. Note that the renormalization of 
(fp /v^ is of higher order in the electromagnetic coupling constant. 

The renormalization of Cj* and Oi is discussed in Ref. [T7] and is given by 

Ci = Zr/C?, Oi = ZnO?, Z^,= (l--^\ \ (2) 

V as e J 
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where /3(as) is the /3-function of standard {ni = 5)-flavor QCD: 

Ai^X^— = , (3) 

with /3(as) = — (as/vr)^/3o + ... and /3o = 11/4 — n;/6. Note that here and in what follows, 
ttg denotes the strong coupling constant in standard five-flavor QCD; it is a function of the 
renormalization scale fi: 

as^a'i'\^,). (4) 



According to Eq. ([T]), the calculation of the gluonic Higgs production cross section or the decay 
rate splits into two parts: On the one hand, the coefficient function Ci has to be determined to 
the desired order. This part is independent of the process considered; the result will be useful 
for any process to which the effective Lagrangian approach is applicable. The second part of the 
calculation is to evaluate the production or decay rate on the basis of the operator Oi within 
five-flavour QCD. 

In the SM, the coefficient function Ci is known to order . The ag-terms have been computed 
in Refs. [18l [7], while the contributions of order were obtained in Refs. [19l [T7] (see also 
Ref. [20]), and the terms of order in Ref. [T7]. For later reference, let us quote the NLO result 
at this point: 

Within the SM framework, there exists an all-order low-energy theorem |17j which expresses 
Cf^ in terms of (3 (as) and 7m(as), the anomalous dimensions of ag and Mt, respectively. The 
result of Ref. [15] for the squark-contribution to the effective theory at order was derived 
by replacing the /3-function with its pure squark contribution, and with the anomalous 
dimension of the squark mass. The obtained expression was multiplied by a factor oc {Mt/Mtf', 
unaffected by renormalization^ 

Our explicit calculation shows that for the correct effective Higgs-gluon coupling, the top squark 
contribution can not be considered on its own. In fact, as soon as the proper renormalization 
of the Higgs-squark coupling is taken into account, the pure squark contribution receives ultra- 
violet divergences which are canceled by a contribution coming from top quarks and gluinos, 
like the ones in Fig.[2](c), [d). As a consequence of this important conceptual difference, also the 
numerical results for hadronic Higgs production have to be re-analyzed. 

In this paper, we follow the diagrammatic approach of Refs. [191 CZ] in order to evaluate the 
effective Higgs-gluon interaction in the MSSM, including terms of order a^. As explained above, 
this involves three types of two-loop Feynman diagrams: Gluonic corrections to top quark 
loops, e.g. Fig.[2](a); gluonic corrections to top squark loops, e.g. Fig.[2](5); gluino-contributions, 
e.g. Fig.[2](c), (d). The resulting effective Lagrangian can then be used to evaluate the Higgs 
production rate at a hadron collider to NLO. 

^Assuming this (non-supersymmetric) theory of scalar quarks, we reproduced the result of Ref. [T5] in our 
diagrammatic approach. Only gluonic corrections to Fig.[l](6) contribute in this case, e.g. Fig.[2](6). 
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(a) (6) 
Figure 1: One- loop Feynman diagrams contributing to gluon fusion in the MSSM. 




(a) 



t , 



(c) 



Figure 2: Sample two- loop Feynman diagrams contributing to gluon fusion in the 
MSSM. Apart from the pure QCD corrections to the top quark (a) and top squark 
(b) triangle, there are also corrections involving the gluino (c), (d). 



2.2 Calculational details 

It is well known that Dimensional Regularization (DREG) |21j breaks supersymmetry. This can 
be cured by adding finite counterterms that restore the super symmetric Ward identities. Alter- 
natively, one can choose a regularization procedure which preserves supersymmetry. A viable 
method which is very similar to DREG is Dimensional Reduction (DRED) [22] (see also [23j). The 
essential difference to DREG is that the tensor algebra is performed in four instead of D = 4 — 2e 
dimensions, whereas the loop integrals are still D dimensional. We checked that the final result 
is the same, both in DRED and in DREG, provided the aforementioned finite counterterms are 
taken into account. 

The generation and evaluation of the Feynman diagrams proceeds fully automatically through 
the chain 

QGRAF Q2E -S- EXP MATAD/MINCER 

m m m m/m 

, , . asymptotic . , , . 

generation analyzation . calculation 

expansion 

Q2E translates the output of QGRAF into a format that is suitable for further manipulation. EXP 
applies asymptotic expansions on a diagrammatical level, according to a user-defined hierarchy of 
mass (possibly also momentum) scales, and produces output that can be given to the integration 
packages MINCER and MATAD immediately. The latter use the Integration-hy- Parts algorithm [29] 
in order to evaluate the single-scale integrals analytically. 



4 



2.3 The effective Lagrangian at next-to-leading order 



Various methods for the evaluation of Ci are described in Refs. \n\ I30j . We fohow the most 
direct one, which means to apply the projector 

to the ggcf) vertex diagrams (sample diagrams are shown in Figs. [T] and [2]), and evaluating them 
at pi = P2 = 0. Here, /i, v and a, b are the Lorentz and color indices of the external gluons, and 
pi , p2 are their momenta. This will result in the product C3 C*? > where is the bare decoupling 
constant of the gluon field obtained from the one-loop gluon propagator involving top quarks, 
top squarks, and gluinosH 



Cl = l + ^{j^+L{e)]+0{al), with 



(7) 



and C2 = vr^/G. Ml/r is the mass for the supersymmetric partner of the left/right-handed top 
quark, Mg is the gluino mass, and Mf the top quark mass. One-loop renormalization of the 
strong coupling constant and the masses (see, e.g., Ref. [H]) leads to Ci[as), where dig is the 
renormalized coupling constant in the full theory. It is related to a^, the coupling in standard 
five-flavor QCD (see Eq. by 

as = {Cgfas, where Cl = - —L{e) + 0{al) . (9) 

L(e) is defined in Eq. ([7|. 

For brevity, we will focus on the following two cases in this letter: 



{A) Mg:^ ML = MR = Mt^ Mt, 
(B) Mr :^Mg:> Ml ~ Mt , 



(10) 



These cases are motivated by the fact that the supersymmetric (SUSY) contributions to the 
gluon-Higgs coupling are proportional to Mt^ / Mt^. Sizable effects are thus only expected if at 
least one of the top squarks has a mass of the order of the top quark mass. Bottom squark 
contributions are suppressed by (Mft/M^)^ • tan/3 and can savely be neglected for not-too-large 
values of tan/3, as we assume them in this paper. Furthermore, as already mentioned in the 
Introduction, we do not consider the mixing of the top squarks. 

The one-loop calculation of Ci involves only diagrams with a single mass scale: Mt or Mt. 
At two-loop order, up to three different masses can appear in the loop diagrams (see, e.g., 
Fig.[2](c), (d)), which makes the calculation quite tedious. However, adopting the hierarchies of 
the cases {A) and {B) introduced above, one can apply asymptotic expansions [32| (see also 



We refrain from quoting terms proportional to 7E and ln47r that drop out of DR-renormalized quantities. 
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Ref. [33] for pedagogic examples) in order to reduce the multi-scale to single-scale Feynman 
diagrams. This leads to significantly simpler integrals and final results of a handy structure 
(powers and logarithms of the masses). 

For later convenience, we parametrize the SUSY contributions to the coefficient function in the 
following way: 



l + ^C^USY 



vr 



where 
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Co{xtL,XtR) = 1 + ^ i^tL + Xtn) , 



XtL 



Ml 
Ml 



XtR 



Ml 
Ml 



(12) 



represents the coupling of the Higgs boson to top quarks which is specified below. We evaluate 
case {A) of Eq. (llOp by first considering the following subcase: 



Mg:$> Ml = Mr = Mt:> Mt. 



(13) 



Expressing the top quark and squark mass in the on-shell scheme, one gets: 



^SUSY 



[A'] 



Xts 



Co{xts,Xts) 



+ xfg a*-^-* 



+ 0{x%) 



(14) 



with 



.(2) 



11 



1 



^ 4 



157 
""36" 

+ Xts 



12 



-L 



1 

18 
1 



11 /14 4 

+ —Ltg + xts I — + -L 



12 



tg 



gs^tg 



37 



+ C2 



12 



12 



-L 



1 



tg 



-,LgsLtg 



V 9 



tg 



25 
12 



+ 3L 



tg 



(15) 



Xts 



Ml 



Xtg 



Ml 

Ml 



L 



tg 



In-^ 

M2 



L 



gs 



M2 
In^ 



(16) 



In the numerical analysis of Section [3l the coefficients a^*) are included up to i = 6. Note that 
the explicit /i dependence drops out which has to be the case as the anomalous dimension of Ci 
has to cancel the // dependence of as in the one- loop term (cf. Eq. ([2])). This is a strong check 
for our calculation. 

An observation concerning Eq. (jl5p is that the two-loop contribution of order Mt'^/M- depends 
logarithmically on Mg which is a consequence of the simultaneous decoupling of all heavy par- 
ticlesH However, in the limit that all SUSY particles are heavy, as expected, the SM result is 
recovered. Let us add that the bare two-loop result even has contributions oc M^^ . Mt^ /Mt^, 
but those are canceled by the on-shell counter term of the squark mass. 

Another remarkable observation concerning Eq. p5p is that higher order terms in Mi' jm^ iden- 
tically vanish, i.e., the coefficients a^'\ i = 0, 1, 2, . . . are exact. This means that Eq. ()15p is valid 



^Note that the Appelquist-Carazzone decoupling theorem [34] is not apphcable here since removing the gluino 
from the MSSM leads to a non-renormalizable theory, if the supersymmetric relation between the top and the 
stop Yukawa coupling should be preserved. 
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for arbitrary values of Mf and Mt, provided that Mg is sufficiently large. In particular, it covers 
case (A) of Eq. ([ini): 



„SUSY 



(17) 



This equality has been checked up to order and it is suggestive that it holds to all orders 

in the inverse gluino mass. 

In order to quantify the condition Mg ^ (Mt, Mt), we show in Fig. [3] the expansion of Eq. (jlSp . 
including successively higher orders in Xgg = Mf/Mg, for fixed values of xts- One clearly 
observes that the expansion breaks down dramatically beyond a certain value of Xgg- However, 
we conclude that for Mg > 2.5 Mt and Mt > Mt, the expansion should be accurate to better 
than 1%. Note also that, as mentined above, the expression diverges logarithmically like ln{xsg) 



as X 



0. 



As an interesting special case, let us consider the limit M^ = Mr = Mt which leads to the 
compact formula 



Ci[ML = MR = Mt, 



-9t 



2^ 



1 + 



a. 



vr 



11 4 



(18) 



It is worth noting that the coefficient a^^^ in Eq. (jlSp is positive only for Mg < e^^^^'^Mt ~ 2.8Mt. 

Considering case (B) of Eq. (fTOl) . we proceed in complete analogy to case (A), i.e., we first 
evaluate the subcase 



{B') Mr > Mg Ml:^ Mt. 
Keeping only the leading term in 1/M/j at NLO, we find 



(19) 



where 



^SUSY 



[B'\ 



XtL 



5(0) 

5(1) 

5(2) 



Co{xtL,XtR 

29 LgR LgL L 



5^ + XLg 6^') + Xlg 6(2) + 0{xtR, Xlg 



48 12 
1 



+ 



— + ^ 
6 2' 



5 ^ 11 r ^13 2 

'18 " 24^^^ + 24^*^ + ""^^^ 18 + 3 ' 
163 1 1 7 ^ 1 ^ 

+ 2^' " 8 " 24^*^ " 4^^^^*^ 



+ XtL 



85 1 3 7 1 r r 1 

36 2 S^'^ ~ 24 ~ 1^9L^i9 ] + ^tL 



+ 0{as) 



49 3 
48 + 2^*^ 



(20) 



(21) 



Ml 
Ml 



Ltg = In 



Ml 

Ml 



XtR 



Ml 

Ml 



XLg 



Ml 

M|' 



Ml 



Ml 



Ml 

M2 



R 



(22) 



We again evaluated the coefficients up to order x^j^ and x^g and use them in the numerical 
analyses below. Note that this result still has a logarithmic dependence on Mr; this comes 
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Mstop = Mtop 



Mstop = 2 Mtop 
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Figure 3: c^^SY q^^^ (-^^ q£ ^j^j^ j^^ ^ ^^^^ (-^^ ^ 2Mt, as a 

function of the ratio Alt / Mg . The individual lines correspond to different orders of 
the expansion in the square of this ratio. The dotted line includes only the leading 
term, whereas the solid line incorporates corrections up to order {Adt/ Mg^^ . 



in addition to the logarithmic dependence on Mg as Mg — )• oo, as it was already observed in 
case (A). 

Also in analogy to case {A), the coefficients ft^*-*, i = 0, 1, 2, . . . are exact, meaning that Eq. (|2ip 
also covers case (B): 



SUSY 



[B] 



SUSY r T3i 



(23) 



For sufficiently large values of Af/j, an analogous study as in case (A) to quantify the condition 
Mg ^ {Ml, Alt) leads to very similar conclusions, so that we refrain from presenting it here. 



3 Higgs Decay Rate and Production Cross Section 
3.1 Decay rate 

From Eq. ([T|) one can derive a general expression for the inclusive (f> ^ gg decay width0 

T{^^gg)=To{l + 5) , 
Fq represents the complete LO result which is given by (see, e.g., Ref. [2]) 

_ Gpa^M^ 



16^/27r3 



E at (A{r,) + -^A{f,L) + j^A{f,n) 
qe{b,t} \ ^^^L ^^\r 



4M2 



2 ' '^qL — 2 ' '^^R 



M^ ' " M. 

A{T)=T[l + {l-T)f{T 



Mr) 



Ml ' 
1 



f [l-f/(f)] 



(24) 



(25) 



Note that at NLO also the final states ggg and qqg contribute to this quantity. 
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Figure 4: The quantity 2c^^^^ as a function of the squark mass for Mg = 1 TeV and 
Mt = 175 GeV. Notice that the SUSY effects are negative. 



with 



/(r) 



arcsm 



In 



1 + 



l-Vl-T 

Relative to the SM case, the top quark coupUng reads 



r > 1, 
r < 1. 



(26) 



9t 



cos a 
sin /3 ' 



9t 



sm a 
sin/3 ' 



(27) 



where h, H denote the hght and heavy neutral scalar Higgs boson, respectively, and a is the 
mixing angle between weak and mass eigenstates in the Higgs sector. Since our focus is on 
not-too large tan f3, we will neglect the effect of bottom quarks and bottom squarks throughout 
the paper. In this limit, gf, and therefore also tan /?, enters our expression as an overall factor. 
At NLO, the correction term can be written as 



a 



vr 



SM 



+ 2d 



,SUSY 



) 



with 



95 7 /11 1 V 

TLi + ni In — ^ 

4 6 V 2 3 7 M2 



(28) 
(29) 



where n/ = 5 is the number of light quark flavors and c^'^^^ is given in Eqs. (|15p and (j21|) . 

The effect of the SUSY corrections is shown in Fig. U] where the quantity 2c^^^^ is plotted as a 
function of the squark mass for Mg = 1 TeV. For Og = 0.1, we get corrections of the order of 
-5%. 
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3.2 Production cross section 
3.2.1 Partonic Results 

The hadronic cross section a^k = a{hk (p + X) for Higgs production can be written as 



) = ^/ dxi dx2 iPi/h{xi)iPj/k{x2)aij 



X1X2 



^ = , (30) 



where i, j denote any partons inside the hadrons /i, A;, and ^p{x) is a parton density, cjjj is the 
partonic cross section, and s is the hadronic cm. energy 

It is convenient to write this partonic cross section as 



Mi 



S 8M| 

where Tq is defined in Eq. (j25p . s is the partonic center-of-mass energy, and 



We find for the NLO terms (g G {u, d, s, c}]f|: 
A« = A«'S^ + 2cSUSY5(l_^)^ 



(31) 



(32) 



a(1),SM 
99 



A 



(1) 



+ 6C2 5(1 - x) + 12 



11 \ ^ rin(l - x) 

l-x 

6(x2 + l-x)2 11 

ln(x) - y(1 - xY 



a(i) = a(i).sm 

99 



:i + (l-x)2)ln- 



(1-x)^ 



12x(-x + x^ + 2) ln(l - x) 



1 + 2x x^ 

3 



A 



(1),SM 
99 



—d-x? 

27U X) , 



(33) 



where c^usY jg (jg^j^ed in Eq. (j29p and the subscript "+" denotes the standard plus distribution. 
A,^j^(x) is the NLO Standard Model expression [6l[7]. All other partonic subprocesses vanish at 
NLO. Note that the SUSY corrections only modify the coefficient of the (5(1 — x) contribution in 
the gluonic channel. 



3.2.2 Hadronic Results 

Fig. [5] shows the NLO K-factor for hadronic Higgs production in the Standard Model and in 
Super symmetry, for two different choices of squark and gluino masses, corresponding to the 
cases (A) and (B). We notice that, even though the overall normalization of the cross section 
in SUSY can be a multiple of the SM value, the QCD corrections are very similar in both cases, 
differing by less than 10%. However, in contrast to Ref. [15], we find that the SUSY effects are 
negative for large values of the gluino mass. 

^ The renormalization and factorization scales have been identified with and can easily be reconstructed 
from lower order results. 
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Figure 5: K- factors at NLO for Higgs productfon in gluon fusion in the Standard 
Model (dotted) and in Supersymmetry. Dashed line: Mr = Ml = Mt = 175 GeV, 
Mg = 1 TeV — Solid line: Mr = 1 TeV, Ml = Mt = 175 GeV, Mg = 400 GeV. 
The bottom Yukawa coupling has been neglected. We use the LO and NLO parton 
densities of MRST2001 [35]. 



4 Conclusions 

In this paper we considered the NLO super symmetric corrections to the production and decay of 
a Higgs boson. We used the framework of an effective Lagrangian where the heavy particles enter 
the coefficient function, Ci, of the operator describing the gluon-Higgs coupling. The practical 
calculation is based on the evaluation of the gluon-Higgs vertex diagrams using asymptotic 
expansions. Our results are in disagreement with Ref. [15j, which is due to the reasons that have 
been discussed above. The SUSY corrections decrease the effects from pure QCD by less than 
5% in the considered parameter space. A more detailed numerical analysis and the inclusion of 
general squark mixing will be presented elsewhere. 
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